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Abstract 

The dynamics of "vicious", continuously growing civilizations (domains), which 
engage in "war" whenever two domains meet, is investigated. In the war event, 
the smaller domain is annihilated, while the larger domain is reduced in size by 



C^"- ' a fraction e of the casualties of the loser. Here e quantifies the fairness of the 



war, with e = 1 corresponding to a fair war with equal casualties on both side, 



■ and e = corresponding to a completely unfair war where the winner suffers no 



casualties. In the heterogeneous version of the model, evolution begins from a 



specified initial distribution of domains, while in the homogeneous system, there 
is a continuous and spatially uniform input of point domains, in addition to the 
growth and warfare. For the heterogeneous case, the rate equations are derived 
and solved, and comparisons with numerical simulations are made. An exact 
solution is also derived for the case of equal size domains in one dimension. The 
heterogeneous system is found to coarsen, with the typical cluster size growing 
linearly in time t and the number density of domains decreases as 1/t. For 



the homogeneous system, two different long-time behaviors function 
of e. When 1/2 < e < 1 (relatively fair wars), a steady state arises which is 
characterized by egalitarian competition between domains of comparable size. In 
the limiting case of e = 1, rate equations which simultaneously account for the 
distribution of domains and that of the intervening gaps are derived and solved. 
The steady state is characterized by domains whose age is typically much larger 
than their size. When < e < 1/2 (unfair wars), a few "superpowers" ultimately 
dominate. Simulations indicate that this coarsening process is characterized by 
power-law temporal behavior, with non-universal e-dependent exponents. Some 
of these features are captured by a deterministic self-similar model, for which the 
characteristic exponents can be computed easily. The transition point e = e c = 
1/2 is characterized by slower than power-law coarsening. 



P. A. C. S. Numbers: 01.75+m, 02.50.-r, 89.90+n 
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I. INTRODUCTION 



(a) Background 

Coarsening phenomena underlie a wide variety of physical processes, such as phase 
ordering and spinodal decomposition [1-3], growth of breath figures [4,5], spin dynam- 
ics [6], and foams [7,8]. The latter system is especially interesting, as the basic phenomena 
are readily observable in many everyday situations. The microscopic rules which govern 
the dynamics of individual bubbles are simply given in terms of the geometry and yet 
lead to varied and intriguing macroscopic behavior in the long-time limit. In addition 
to continuous growth and shrinking of individual bubbles, there are discontinuous bubble 
"popping" events which lead to rearrangement on a larger scale than individual bubble 
growth. Another attractive feature of such geometric coarsening processes is that they 
naturally suggest idealizations which may be exactly soluble. One such example, which 
can be viewed as a limiting case for breath figure growth, is the coarsening of an array of 
contiguous domains in one dimension [9-17]. 

The mechanisms that govern domain evolution in these types of coarsening processes 
appear to have natural counterparts in social phenomena. For example, the competition 
between cultures has led to a rich historical record [18] in which certain civilizations are 
dominant for long time periods only to suddenly disappear. Conversely, other civilizations 
persist for very long times even though they are relatively small. Motivated by these basic 
historical facts and by a qualitative appreciation for coarsening phenomena, we introduce 
and investigate a simple, yet relatively general model for war between vicious civilizations. 
In our model, domains grow continuously and an encounter between two civilizations leads 
to a war where the smaller combatant is annihilated, and the larger civilization suffers a 
specified number of casualties. This model exhibits a rich variety of dynamical behaviors, 
including the features of extinction and persistence of civilizations. 

(b) The War Model 

In our model, space is populated by a collection of civilizations, each of which is 
represented by a spatial domain of a particular size (or population). Each domain grows 
either continuously or in discrete steps at a constant rate, in which the boundaries move 
at velocity V. Whenever two civilizations of size i and size j meet (with j > i without 
loss of generality), they engage in a war where the smaller domain is annihilated, while 
the larger domain suffers ei casualties, so that its size changes to j — ei (Fig. 1). We 
define the casualties to occur at the battlefront so that the frontier of the winner retreats 
by a distance ei in one dimension. Here e measures the "fairness" of the warfare event; 
e = 1 corresponds to a fair war, in which the winner and loser suffer the same number 
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of casualties, and e = corresponds to a completely unfair war in which the winner 
emerges unscathed. Subsequently, the survivor civilization continues its growth until the 
next war. We are interested in determining the long-time dynamical behavior and the 
spatial distribution of civilizations under these conditions. Although our model is naive 
and drastic (since the loser is annihilated), the essential features of growth and sudden 
diminishment by war are incorporated. Thus suitable generalizations of our model may be 
appropriate for describing quantitative aspects of social history. 

There are two general conditions under which our model leads to interesting asymp- 
totic dynamics. In the heterogeneous warfare process, the system begins with a distribution 
of nascent domains which subsequently undergo growth and warfare. In the homogeneous 
process, there is temporally constant and spatially uniform input of point domains which 
subsequently grow and experience wars whenever two domains meet. From a social per- 
spective, this input could be viewed as arising from the remnants of destroyed civilizations. 

(c) Overview of Dynamical Properties 

For the heterogeneous warfare process, we have primarily concentrated on the case 
e = 1, although qualitatively similar results are anticipated for any value of e. This system 
exhibits coarsening, in which the domain size distribution approaches a scaling form in 
x/t, where x is the domain size and t is the time. Consequently the average domain size, 
(x(t)), grows linearly in time. These qualitative results hold both in the mean-field limit 
and in one dimension. From the rate equations, the scaling function for the domain size 
distribution is one-half of a sinusoid. In contrast, simulations in one dimension reveal a 
scaling function for the size distribution which is triangular in shape. 

The homogeneous warfare process exhibits a richer phenomenology which is controlled 
by the fairness parameter e. For 1/2 < e < 1, a steady-state arises, which is characterized 
by egalitarian competition between domains of similar size, while for e < 1/2, coarsening 
occurs in which a few superpowers ultimately dominate. The existence of a threshold 
between these two regimes can be qualitatively justified by considering the outcome of a 
war between two neighboring domains with different birth times (or equivalently, sizes). 
For e c = 1/2, when two neighboring domains meet and engage in war, the frontier of the 
survivor after the war will retreat to exactly the same position that existed at the birth- 
time of the smaller civilization (Fig. 2). Thus for e ^ 1/2, the frontier of the victor after 
the warfare event should either advance or retreat, respectively, compared to its position 
at the birth-time of the weaker combatant. This difference suggests the existence of the 
aforementioned steady state and coarsening regimes which are separated by e c = 1/2. 

An interesting feature of the steady state regime, e > e c , is that a non-trivial joint 
distribution of civilization age and size arises. The connection between these two attributes 
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is subtle, since a long-lived civilization, whose age is much greater than its size (in appro- 
priately scaled units where V = 1/2), can arise by a sequence of fortuitous events. To 
characterize this age-size distribution, we have formulated a version of the rate equations 
which simultaneously accounts for the distribution of domains, as well as the distribu- 
tion of gaps between domains. This approach appears to provide some exact results for 
the age-size distribution in one dimension, at least in the analytically tractable case of 
e = 1. For the distribution of domain sizes (integrated over all domain ages), n(x), we find 
n(x) = e~ x / x " with xq = l/^fjl, in excellent agreement with simulations. Here \i is the 
rate per unit length at which new civilizations are introduced. More interestingly, the age 
distribution (integrated over all domain sizes), n(r), has an exponential tail, n(r) oc e _r / T ° 
as t — > oo, but with an anomalously large characteristic age tq which is much larger than 
the typical size. Simulations in one dimension give an even larger value for the ratio tq/xq, 
a result which can be attributed to the existence of significant anti-correlations in the ages 
of neighboring domains. For the case of general e above the threshold, e c < e < 1, our 
simulations indicate that basic quantities, such as the average domain size, domain age, 
and coverage, all eventually reach steady-state values. Correspondingly, both the domain 
size and age distributions decay exponentially in this regime. 

In the complementary regime of e < e c , warfare events sufficiently favor the victor that 
superpowers ultimately emerge. The existence of these large domains strongly modifies the 
effect of the continuous input of small domains, so that a steady state is not achieved. Our 
simulations indicate a continuous coarsening of domains in which basic time-dependent 
quantities, such as the mean domain size and the number density of domain, exhibit non- 
universal e-dependent power law behavior in time. In particular, for the extreme limit of 
e = 0, where the victor in a war suffers no casualties, the number of domains and the 
density of empty space decay as t~ J in one dimension, with 7 « 1/3. A deterministic 
idealization of the war model is introduced which provides a useful description of this 
coarsening phenomenon. As e approaches e c from below, coarsening becomes slower, as 
evidenced by the e dependence of the exponents. For e ~ e c , marginal behavior occurs in 
which the average domain size grows extremely slowly, while the average age continues to 
grow as a power law in time. 

In Sec. II, we study the heterogeneous war model with fair wars, e = 1. Results from a 
mean-field theory and numerical simulations are presented. We also treat the special case 
where domains all have the same initial size, so that both domains disappear in a warfare 
event. This limiting situation turns out to be exactly soluble in one dimension for both 
the asymptotic value of the coverage, as well as the domain size distribution. We next 
turn to the homogeneous process in Sec. III. From the rate equations, the equation of the 
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steady-state solution for the joint age-size distribution is obtained, again for the particular 
case of fair war, e = 1. The extension of the rate equations for e < 1 is also discussed and 
an argument for the existence of a transition in the kinetic behavior at e = e c is presented. 
We next present simulation results in one dimension, focusing on general properties as 
a function of the fairness parameter e. Basic features of the e > e c steady-state regime 
and the e < e c coarsening regime are outlined. A deterministic self-similar model is then 
introduced to help understand the coarsening process for e = 0. In Sec. IV, we give a brief 
summary and discuss several extensions of our model. Various calculational details are 
given in the Appendices. 

II. HETEROGENEOUS WAR 
(a) Mean-Field Theory 

Consider the mean-field limit, in which pairs of domains are randomly picked to un- 
dergo warfare. With this interaction rule, the evolution of the domain size distribution is 
described by the rate equations [19], 

CO CO 

Ck(t) = ^2ci(t)c l+k (t) - c k (t)^2ci(t) + A(c fc _i(t) - c fc (t)). (1) 
i=l i=l 

Here Ck(t) is the concentration of domains of size k at time t. The first term accounts for 

the gain of /c-domains due to a war between domains of size i and i + k. Similarly, the 

second term accounts for the loss of /c-domains because of war between a domain of size 

k and any other domain. The gain and loss of /c-domains due to constant growth at rate 

A are described by the last two terms. We have implicitly assumed that the warfare rate 

is independent of the combatant sizes. Eq. (1) also assumes a minimal size of unity, so 

that all sizes are integers. The extension to the continuum case is straightforward; and 

this description will be employed for the homogeneous war model. 

For the system described by Eq. (1), the average civilization size grows linearly with 

time and the size distribution is one-half of a sinusoid. To derive these results, we first 

identify the appropriate scaling variable. For this purpose, let us temporarily neglect the 

effect of war. The size distribution for such peacefully flourishing civilizations is found 

from the rate equation tk = X(ck-i — Ck). This is just the Poisson process, with solution, 

for the "Adam" initial condition of Ck(t = 0) = 8k,i, 



<*(«> = fcV ~ w 



(k - Xtf 



While the exact result is specific for the monodisperse initial condition, the Gaussian 
approximation provides universal asymptotics for an arbitrary compact initial size distri- 
bution. Thus for peaceful civilizations, the size distribution is peaked around k = At, with 
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dispersion yXt. Clearly, warfare produces sizes with k < Xt, so that the size distribution 
for warring civilizations should be no n- vanishing for all k < Xt. This suggests that the 
appropriate scaling variable is x = k/Xt, with < x < 1. The leading region where k ~ Xt 
can be expected to have a fine structure over an extent of the order of VXt, as in the 
special case of peaceful civilizations. In the following, we ignore this detailed structure. 

To solve the rate equations, it proves useful to consider first the civilization number 
density, N(t) = J2k>i c k(t), which, from Eq. (2), satisfies 

1 00 1 

N(t) = --J2c k (tf--N(tf. (3) 

k=i 

Asymptotically, the first term on the right-hand side is clearly negligible; therefore, as 
t — > oo, N(t) ~ 2/t. We thus expect that scaling form for the civilization size distribution 
is 

1 k 

c k (t) = ^ C(x), with x = — . (4) 

The time dependent prefactor in Eq. (4) guarantees that N(t) oc t~ x . Furthermore, the 
relation N(t) ~ 2/t is quantitatively satisfied if dxC(x) = 2. 

Substituting the scaling ansatz into the rate equations, one finds that in the continuous 
limit, the scaling function satisfies an integro-differential equation of an anti-convolution 
form, 

(l-x)C'(x) = [ dyC(y)C(x + y), (5) 
Jo 

where the prime denotes the differentiating with respect to x. We are unable to solve this 
equation by a systematic approach, and therefore resort to trial- and-error. Since C' > 0, 
C(x) is monotonically increasing in x. On physical grounds, we expect that C(0) = 0, while 
analysis of Eq. (5) in the vicinity of x = 1 gives 

C(x)=C(l)-(l-x) 2 ^M + ..., (6) 

i. e., C'(l) = and C"(l) < 0. Polynomial test functions fail to satisfy Eq. (5). (For a 
polynomial, say, of degree g 1 the left-hand side of Eq. (5) is a polynomial of degree g while 
the right-hand side has degree 2(7 + 1). We find, however, that the simplest appropriate 
transcendental test function, C(x) = n sin( 2 rf), solves Eq. (5). (The constants in this 
expression for C(x) are chosen to satisfy C'(l) = and the sum rule for C(x).) Combining 
Eq. (6) with the scaling form of Eq. (4) gives the basic result, 

^sin(^), k<Xt, 
c k (t) * < (7) 
0, k > Xt. 
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From this solution, the moments of the size distribution, M n (t) = J2k>i k n c/c(£) are 
given by 

2 4A 8(tt-2) „ 

M --, M 1 = — , M 2 = v ; A 2 t, . . . , 

/on+l \ ( 8 ) 

M n =f—^y da;x n siiixj A n t n_1 . 

Thus, e. 5., M (£) is the civilization number density, M 1 (t) is the mass density, or coverage, 
Mi(t) = M(t); etc.. From Eq. (8), any reasonable measure of the typical domain size, e. g., 
[M n (t)/M (t)] 1 / n , increases linearly with time. 

(b) Simulation Results in One Dimension 

A straightforward way to simulate our warfare model in one dimension is by molecular 
dynamics. Within a continuum and deterministic description for domain growth, one 
identifies the minimum of the conflict times between all pairs of nearest-neighbor domains 
at any given stage. The system evolves freely until this minimum conflict time, at which 
point a war occurs, where one domain (the smaller) disappears and its (larger) nearest 
neighbor shrinks. After this event, all pairwise conflict times are recomputed and the 
overall update process is repeated. While simple to implement, this molecular dynamics is 
relatively inefficient, as the computation time is proportional to the square of the number 
of domains. We therefore employed an alternative algorithm which leads to a savings of 
almost two orders of magnitude in time for a system with 10 5 initial civilizations, compared 
to molecular dynamics [20]. In our approach, we first determine the conflict times for all 
neighboring civilizations and sort them in ascending order. Instead of re-computing conflict 
times after each war, we continue to use the pre-sorted times for carrying out successive 
warfare events until a threshold is reached. This threshold is determined by first computing 
the new nearest-neighbor conflict time that is created as a result of the current war and 
comparing this new time, as well as the nearest-neighbor conflict times which are "lost" by 
the current war, with the next conflict time on the pre-sorted list. If any of these putative 
times are less than this next pre-sorted time, an inconsistency would arise at the next step. 
It is then necessary to re-compute and re-order all conflict times. This exhaustive molecular 
dynamics step needs to be performed relatively rarely, leading to considerable saving in 
simulation time. It is possible, at the expense of algorithmic simplicity, to eliminate the 
molecular dynamics step entirely by constantly ordering the list of conflict times as conflict 
times are created and destroyed in each warfare event. 

For a polydisperse initial distribution of civilization sizes, our simulations show that 
the time-dependent size distribution evolves to a nearly universal scaling form in the asymp- 
totic limit. Details of the initial condition are irrelevant as long as they are not singular in 
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character. Our numerical results are based on using a Poisson initial distribution for both 
the sizes of civilizations and the intervening gaps. However, the shape of the size distri- 
bution is influenced by details of the warfare event, such as the location of the removed 
portion of the victorious domain. In our simulations, this removed portion is adjacent to 
the battlefront (Fig. 1). Other rules are possible and perhaps natural, e. g., one could 
define a rule in which the center-of-mass of the survivor remains fixed after the warfare 
event. Because of this detail dependence of the size distribution, its quantitative char- 
acterization is of limited value. For the casualty rule adopted here, the distribution has 
simple triangular shape (Fig. 3). 

An interesting feature from the simulations is that the sizes of nearby domains are 
virtually uncorrelated. That is, the size correlation function, C s (r) = (siSi+ r )/ (si) 2 — l ~ 
for r > 1, where Si is the size of the i th domain. For r = 1, the simulations give C s (l) ~ 
—0.01 which is at least 3 times larger than the correlation function for any other value of 
r. Thus two large civilizations are less likely to coexist peacefully as nearest neighbors; 
rather, a large domain is slightly more likely to be surrounded by small neighbors and vice 
versa. 



For heterogeneous war in one dimension, the case of initial equal-size civilizations 
is unique because both combatants are eliminated in a war and the equal-size distribu- 
tion is preserved. This case turns out to be exactly soluble by appealing to a connection 
with domain coarsening processes (see, e. g., Refs. 9-13), and generalizing the approaches 
in Refs. 14-17. In domain coarsening, which we may view as being "dual" to war for 
equal size domains (Fig. 4), the system consists of contiguous domains of arbitrary sizes 
and coarsening occurs by successive elimination of the smallest domain. When the walls 
associated with this minimal size domain disappear, the other walls remain fixed. This 
successive domain elimination corresponds exactly to the pairwise annihilation of the two 
closest domains in the war model, as illustrated in Fig. 4. Note also that since all civiliza- 
tions have the same size, their growth rate is immaterial, and models with size-dependent 
growth may be solved by the same approach as that used for size independent growth by 
using the domain length L as the time parameter. 

Let n(l, L) dl be the number of neighboring civilizations of size L whose centers are 
separated by a distance which is within [/, I + dl]. Using the equivalence to coarsening, we 
term the interval between the centers of neighboring civilizations a "domain" . The total 
number of such surviving domains is 



(c) Equal Size Domains 




(9) 
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It proves useful to normalize this quantity, f(l,L) = n{l,L)/N{L), and then define the 
(almost) scaling form, F(x,L) = Lf(l,L), with x = l/L. The absence of correlations [15- 
17] between domains in the dual coarsening process is crucial since it implies that mean-field 
rate equation for F(x, L) is exact. This rate equation reads [14] 

L ^F(x, L) = F{x, L)+x-^-F(x, L) + 6(x-S)G(L) J dyF(y, L)F(x-y-l, L), (10) 

where G(L) = F(x = 1,L). Eq. (10) can be obtained by a straightforward enumeration 
of the outcomes that arise from the elimination of the smallest domains (see, e. g., the 
derivation of Eq. (9) in Ref. 17). For example, the last term in the right-hand side of 
Eq. (10) describes the formation rate of an x-domain by elimination of the smallest domain, 
of scaled length 1, which is situated between two domains of scaled lengths y and x — y — 1. 
The step function 9(x — 3) ensures that the resulting domain will be at least 3 times larger 
than the minimal domain. 

From the solution to Eq. (10) in the long-time limit, or equivalently, L — > oo, the 
asymptotic coverage is 

M 00 = ^^ 0.28073, (11a) 

where 7 = 0.5772156 is Euler-Masceroni constant. (This asymptotic solution, first given 
in Ref. 14, as well as the full time-dependent solution are detailed in Appendix A). Addi- 
tionally, the number density N(L) of these equal-size civilizations asymptotically is 

N(L) ~ MooZT 1 . (12) 

These behaviors for the number density and the coverage are qualitatively similar to the 
corresponding mean-field results. Additionally, from the complete time-dependent solu- 
tion of the rate equations in Appendix A, one also finds the asymptotic expansion of the 
coverage, 

M(i) =2? + I + '- <116) 
i. e., the first dominant correction decays as L~ l . Here L is a measure of the physical 
time, and the coefficient A depends on the details of the initial size distribution. 

While it would be interesting to investigate the domain evolution in heterogeneous was 
for general < e < 1, the equal-size property is lost as the process develops, and an exact 
treatment does not seem possible. However, the extreme case of completely unfair war, 
e = 0, still enjoys the property that an equal-size distribution remains invariant during 
the evolution if one defines that one of the combatants (picked randomly) is annihilated 
in a war while the other remains unchanged. Such a model may be solvable by techniques 
similar to those employed in the above case of equal size domains engaging in fair war. 
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III. HOMOGENEOUS WARFARE 

We now consider the effect of a temporally and spatially homogeneous input of size- 
less civilizations on the dynamics. As discussed in the Introduction, two fundamentally 
different long time behaviors can occur, depending on the value of the fairness parameter 
e. For < e < 1/2, power-law coarsening occurs, leading to the emergence of a few super- 
powers. Conversely, for 1/2 < e < 1, a steady state arises, with egalitarian competition 
between comparable size domains. For the latter situation, it is plausible that a mean- 
field approach might be accurate, since the input leads to a well-mixed state. The rate 
equations of the previous section are not suitable, however, since the restriction to nearest- 
neighbor interactions is not accounted for. Our goal here is to construct rate equations for 
the driven one-dimensional process which incorporates the obvious restrictions associated 
with one spatial dimension. We are able to solve for the steady state of these governing 
equations in the fair war case of e = 1. 

(a) Rate Equations for Fair Wars 

In one dimension, civilizations are represented by non-overlapping intervals, with the 
civilization size equal to the interval length. It is now convenient to assume continuous 
and deterministic civilization growth in which boundaries move with constant velocity V. 
Thus the random birth times and placement of new civilizations (whose initial size may 
be taken as zero) are the only sources of randomness. 

Consider the fair war case of e = 1. To write the rate equations, we first introduce 
the distribution functions, n(x,t) and m(x,t), which are, respectively, the density of civi- 
lizations of size x and the density of inter-civilization gaps of size x at time t. The number 
density of civilizations can be written equivalently as 

POO POO 

N(t) = / dxn(x,t) = / dxm(x,t). (13) 
Jo Jo 

Thus the fraction of covered space is M{t) = J °° dx xn{x,t), while the fraction of empty 
space is E(t) = J °° dxxm(x,t), with M(t) + E(t) = 1. 
The rate equations for n(x,t) and m(x,t) are 



d n(y,t) n(x + y,t) n(x,t) 



o 



f d d\ f°° 

I — - 2 V— J m(x, t) = 2ji J dy m{y, t) - iixm(x, t) 



N(t) N(t) N(t) 

2Vmo(t) m{x, t) 



+ fjL6(x)E(t), 

(14a) 



+ 



N(t) 

(146) 

/ dzm(z,t)n((x — z)/2,t) / dyn(y,t). 

Jo J(x-z)/2 



N 3 {t) 
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In these equations, m (t) = m(x = 0, £) is the density of gaps of size zero and \i is the 
birth rate of new domains per unit length. The spatial derivative term in these equations 
accounts for the continuous growth of civilizations, in Eq. (14a), and the shrinking of gaps, 
in Eq. (14b). The right-hand sides account for the evolution as a result of interactions 
(Fig. 5). The first term on the right-hand side of Eq. (14a) gives the production rate for 
domains of length x as a result of a war between domains of size y and x + y. Such an event 
occurs only when the gap between these two domains vanishes - hence the factor mo(t). 
In the mean-field approximation, the rate for this process is proportional to the product of 
mo(t) and the probabilities n(y,t)/N(t) and n(x + y,t)/N(t); the factor 4V accounts for 
the two possible locations of the combatants, (y, x + y) and (x + y,y), times the rate 2V at 
which the gap vanishes. The second term on the right-hand side of Eq. (14a) accounts for 
wars between an x-domain and an arbitrary size right- or left-neighbor. The last term gives 
the rate at which size-less civilizations are created in empty space. Only this last term is 
exact a priori, because there is no factorization of multi-particle correlation functions into 
single particle densities. 

The terms on the right-hand side of Eq. (14b) are explained similarly. The first 
two terms arise from the "fragmentation" of an empty interval due to the input of new 
civilizations. The gain term accounts for the production of an x-gap due to the two ways 
in which a y-gap can be broken into a gap size y — x and x by the input. The second term 
accounts for the loss of x-gaps due to their total rate of breakup as a result of the input. 
These two terms are, again, presumably exact. The last two terms describe how m(x,t) 
evolves by war. The loss term arises because a war, which is adjacent to an x-gap, leads to 
the removal of the x-gap if the adjoining civilization is the loser. There is a cancelation of 
a factor of 2, to account for the two possible locations of the warfare event, with a factor of 
1/2, to account for the possibility that the loser may not be adjacent to the gap. The total 
rate for any war, independent of the size of the combatants, is simply 2Vmo(t). Finally, 
the gain term arises from wars between a civilization of size (x — z)/2 and y > (x — z)/2 
which is adjacent to a z-gap. Since there are x — z casualties in the war, the initial gap of 
size z grows to size x. 

(b) Steady State Properties for Fair Wars 

To determine the steady state properties of these equations, it is first helpful to con- 
sider the rate equation for the total number density of civilizations, 

^M = -2Vmo(t) + fiE(t). (15) 

This equation is exact; it can be derived directly on physical grounds and also follows by 
integrating Eq. (14(a)) over all x. In the steady state, Eq. (15) becomes 2Vmo = \iE. 
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Substituting this into the steady-state version of Eq. (14(a)) gives 

n(y) n(x + y) n(x) 



d , . 

— nix) = IrriQ 
dx 



Jo 



dy 



N N N 



+ m S(x). (16) 



We seek an exponential solution to this equation. The presence of the 5-function 
implies that n(x = 0) = m®. Thus we hypothesize that n(x) = m§e~ ax . It is easily 
verified that this satisfies Eq. (16). Additionally, from 2Vmo — jxE = fj, J °° xm(x) dx and 
n(x) = moe~ ax , the parameter a is determined by 2Vmo = [i (l — ^r)- To complete the 
solution, we now consider the rate equation for m(x). It is again natural to attempt the 
same exponential form for the gap distribution, m(x) = moe~ ax . Substituting this ansatz 
into the steady version of Eq. (14(b)), consistency is achieved if \i = 2Va 2 . Combining 
with the previous relation gives mo = fJ>/4V. Thus we finally arrive at the steady state 
solution 

n(x) = m(x) = exp (^ x \Jl^j ■ ( 17 ) 

This gives the steady state number density, iV = ^/j,/8V and coverage, M = E = 1/2. 
Results from our numerical simulations of the homogeneous warfare model in one dimension 
are indistinguishable from these results, suggesting that this rate equation approach gives 
the exact the size distribution. 

Consider now the civilization age distribution. This turns out to be a considerably 
more interesting but subtle characteristic of the steady state. Although the age distribution 
is asymptotically an exponentially decaying function of age, the characteristic age is much 
larger than the naive expectation of the characteristic size divided by V = 1/2. Thus a 
typical civilization survives many wars before it is ultimately extinguished. To determine 
the age distribution, it is helpful to consider the more fundamental steady-state joint age- 
size distribution, n(x, r), defined as the density of civilizations of size x and age r, (with 
x < 2Vt). From this joint distribution, the steady-state size distribution of civilizations 
of any age is clearly given by 

/>oo 

n(x) = / drn(x,r), (18a) 

Jx/2V 

while the steady-state age distribution of any size civilizations is given by 

p2Vr 

n ( T ) = / dxn{x,T). (186) 
Jo 

By the nature of the warfare process, the joint age-size distribution consists of two 
components, 

n{x,r) = M(x,t)+1(t)5(x-2Vt). (19) 
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The first term accounts for "mature" civilizations which have experienced at least one war. 
The size of such civilizations is strictly less than the maximum possible size, x max (r) = 
2W, at a given age r. The second (singular) term accounts for "innocent" civilizations 
which have not experienced any war during their lifetimes. The 5-function ensures that 
these innocents are at the maximum size for a given age. These two components of the age- 
size distribution obey different rate equations. The equation for the density of innocents 
is readily soluble and this facilitates the full solution. 

The rate equation for the density of innocent civilizations is 

^ = -4V^l(r) + ^E5(r). (20a) 

The two terms account for the net change of innocent civilizations by warfare and input, 
respectively. Here the fraction of empty space in the steady state, E = 1/2, as derived 
above. Solving Eq. (20a) yields X(r) = |^exp(— T\fd>Vjj), i. e., the age distribution of 
innocent civilizations is purely exponential. 

The rate equation for the density of mature civilizations is 



— + 2V—)M(x,T)=4Vm 



~ x/2V , n(y) M(x + y,r) M(x,t 
dy 



N N N 

J(r) n(2Vr - x) 



(206) 



+4Vm ~- 



N N 



where N = ^/j,/8V is the steady state civilization density. The first term on the right- 
hand side accounts for the gain of mature civilizations size x and age r due to a war 
between a mature (x + y, r)-civilization and one (either mature or innocent) of size y and 
arbitrary age. Similarly, the second term accounts for the loss of an (x, r)-civilization due 
to its undergoing warfare. The last term accounts for the creation of (x, r) -civilizations 
due to an innocent civilization of size x and age r = x/2V experiencing its first war with 
a civilization of size 2Vt — x and arbitrary age. 

The full steady state civilization density n(x) and the density of innocent civilizations 
X(r) are already known. Therefore, Eq. (20b) is a linear integro-differential equation with 
non-constant coefficients. It is possible to reduce Eq. (20b) to the Klein-Gordon equation. 
This reduction, as well as the solution to the resulting boundary value problem, is detailed 
in Appendix B. From this solution, we find the steady-state age distribution 

n(r) ocr- 3/2 exp[-(3- v / 8)r v / 2V>] (21) 

in the large age limit. Therefore, the characteristic age is a factor (3 - v^)" 1 « 5.82 
times larger than the characteristic size divided by the rate of growth 2V . Thus a typical 
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domain survives approximately 6 wars before its ultimate death. Also, it is possible to show 
that the average cluster age (r) = J °° JJ" dr dx r n(x , r) = 2/^/2VjI. On the other hand, 
the constant input at unit rate implies that the domain life expectancy equals \j \J2V [i. 
This inequality between life expectancy and average age is what may be expected in the 
harsh environment defined by our war model. Domains are especially vulnerable close to 
their time of birth, but become progressively more robust as they grow. This behavior 
is akin to that of sea turtles which are most susceptible to predation immediately after 
being hatched. However, if the hatching survives its initial trek into the ocean, it has a 
reasonable chance of living to an old age. 

(c) Steady State Properties for Nearly Fair Wars 

In analogy to the case of exactly fair war, we expect that our equation approach should 
provide an accurate description for steady-state properties in the regime of nearly fair war, 
e < 1. The rate equations in this case are straightforward generalizations of Eqs. (14) 
which, in the steady state, become (taking V = 1/2) 



dn 

— = 2m 

dx 



L 



! *' n(y) n(x + ey) n(x) 



dy 



N N N 



+ m 5(x) 7 (22a) 



— = nxm(x) - 2fj, / dym(y) + m — dz m(z)n f J l_^dyn(y). 

(226) 

While we have been unable to solve these equations, it is possible to show that the 
character of their solution changes as e decreases from 1 to 0. This suggests that steady 
state exists only for a limited range of e > e c , while a different type of solution exists 
otherwise. For this purpose, it is sufficient to consider Eq. (22a) . Using the transformation 

n(x)=moV(0, C=^, (23) 



one can rewrite the rate equation as 

ldv r^ 1 -^ 

2 d£ 



/ d V V( V )V^ + en)-V(0, (24) 

Jo 



which is to be solved subject to the boundary condition P(£ = 0) = 1. For the case of 
completely unfair war, e = 0, we find the explicit solution 

m - ^ (^) 

This solution has a serious flaw in that the first moment of the distribution V(^) 
is divergent. However, the first moment of the unsealed size distribution, n(x,t), is the 
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fraction of covered space and must clearly be finite. Thus a physically acceptable steady- 
state solution does not exist for e = 0. On the other hand, we have previously seen that 
Eq. (24) does admit a reasonable steady solution, P(£) = e~^, when e = 1. We therefore 
conclude that the character of the solution to the rate equation changes for some value 
of e between and 1. Unfortunately, we are unable to determine the threshold value of e 
below which Eq. (24) has no physically acceptable solution. 

It is also worth emphasizing that a rate equation description may not even be appli- 
cable when e < e c . When a steady state exists, there is sufficient empty space available 
in the system for the steady input to act as a relatively effective mixing mechanism. This 
supports the notion that a mean-field rate equation approach could provide an exact de- 
scription of some steady-state properties, as discussed in the following section. On the 
other hand, if the system coarsens, the fraction of empty space vanishes and the input 
becomes progressively less successful in giving birth to new domains. Under this circum- 
stance, it is not evident that the time-dependent rate equations have the potential to fully 
capture the time evolution of the system. 

(d) Simulation Results in One Dimension 

To test our analytical predictions for e = 1 and to map out the dynamical behavior for 
general < e < 1, we have performed molecular dynamics simulations for one-dimensional 
systems with between 1000 to 16000 initial domains, with an input rate of size-less civi- 
lizations per unit length \i = 1/2, and the growth velocity V = 1/2. When a steady state 
arises, any initial condition would be, in principle, suitable. However, to reduce the extent 
of the early time regime, we empirically found that a good choice for the initial condition 
is a Poisson distribution for both the domain and inter-domain gap sizes, each with a 
characteristic length of unity. A summary of our numerical results is given in Tables 1 and 
2. 

For e = 1, our analytical predictions for the steady-state values of the coverage, 
concentration, and also the form of the size distribution are confirmed. In particular, for the 
size distribution, we find n(x) = e~ x / x ° , with xq very close to the exact value of 1/ y/JI = \[2 
(Fig. 6). Simulations for the domain age distribution also suggest that the asymptotic tail is 
exponential, namely n(r) ~ e _T / T °, but with To = 24 ±2, which is to be compared with the 
analytical result from the rate equation r = ((3 — v / 8)y / A*))~ 1 ~ 8.23. The source of this 
discrepancy appears to be the existence of significant correlations in the ages of neighboring 
civilizations, a feature which would render the rate equations inaccurate. This correlation 
arises because the domain age remains unaffected by war, so any age correlations which do 
develop between neighboring domains persist until all of these domains die. On the other 
hand, the domain size is affected by wars, so that correlations in domain sizes should be 
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inhibited by the evolution itself. To check this hypothesis we measured the size and age 
correlation functions, C s {r) = (siSi +r ) / (si) 2 — 1 and C T (r) = (riTi +r ) / (r^) 2 — 1, where Sj 
and Ti refer to the size and age of the i th domain, respectively. As anticipated, C s (r) is 
very close to zero for all r. However, the age correlation function C T (r) is systematically 
negative for r < 5 (Fig. 7), which implies that old civilizations are less likely to coexist close 
to each other. The absence of size correlations and presence of age correlations suggests 
that the mean-field approach should be quantitatively accurate for the size distribution 
but not for the age distribution. 

Although the geometrical properties of domains are time independent in the steady 
state, their properties as a function of domain age are not stationary. This age dependence 
may provide a useful and deeper characterization of the steady state. One such example 
which appears especially intriguing is the behavior of the average domain size as a function 
of the corresponding domain age (Fig. 8). This size grows very slowly and ultimately 
saturates at a finite value as t — > oo which is estimated to be approximately 4. From 
the solution to the rate equations, we find that the average size has the asymptotic form 
(2 + 2^) —const, jr. Thus for a domain to be long-lived, it must not be exceptionally 
large. 

For e c < e < 1, our simulations show that a steady state is eventually reached, but 
that the time needed to attain this steady state grows as e approaches e c from above 
(Fig. 9). Correspondingly, the steady state values of fraction of empty space, E, and the 
concentration of domains, N, become smaller for decreasing e. For example, for e = 1 
and \i = 1/2, our analytical results give E = N = 1/2, while for e = 2/3 and \i = 1/2, 
simulations give E ps 0.35 and N ps 0.25. Additionally, the corresponding domain size 
and age distributions for e c < e < 1 appear to have the same functional forms as in the 
fair-war limit of e = 1. In particular, for e = 2/3 these distributions are n(x) ~ e ~ x / x o^ 
with xq = 5 ± 0.5, and n(r) ~ e _T / r °, with tq = 300 ± 50. This increase in xo and tq as 
e — > e c from above has been expected. As e is decreased, conflicts become less devastating 
for the survivors, so that they may grow larger and live longer. 

When e < e c , the system never reaches a steady state (Fig. 9), and the evolution of a 
finite size system ends when a single superpower occupies the entire space. To quantify this 
coarsening, we consider several basic quantities including the average domain size (x(t)), 
the average maximum domain size (# max (£)), the number density of domains N(t), the 
fraction of empty space E(t), the average domain age (r(t)), and the exponents associated 
with their asymptotic behavior: 

(aw(£)> oc t^ e \ N(t) oc t~^ £ \ E(t) oc t~^ e \ (r{t)) oc . (26a) 
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As written, these power laws are found to be non-universal, with e-dependent ex- 
ponents (Table 1). Further, some of these quantities are interrelated. For example, 
(x(t)) = J dxxn(x,t)/ J dxn(x,t) = (1 - E(t))/N(t). Therefore (x{t)) oc 1/N(t) oc t^ e \ 
Furthermore, the fraction of empty space E(t) and concentration N(t) have the same time 
dependence. To establish this, it is helpful to introduce the normalized domain size dis- 
tribution p(x,t) = n(x,i)/N{t). Numerically, we find that this distribution approaches 
a stationary form, p(x,t) — > p(x), in the long time limit, and with a universal exponent 
(Table 2), 

p(x) oc x~ s when x — > oo. (266) 

This stationarity is, in fact, related to the equivalence between N(t) and E(t). Since 
the concentration of zero size domains is constant and equal to \x for each inter-domain gap, 
and since empty gaps account for a fraction E(t) of the system size, n(x = 0, t) = fiE(t). 
On the other hand, the time dependence in n(x, t) = p(x)N(t) appears only through N(t). 
We therefore conclude that N(t) ~ E(t). It should be noted, however, that this relation 
is valid only if p(x) decreases sufficiently fast for large x, so that J 1 p(x) dx converges as 
t — > oo. Our simulations indicate that this is indeed the case, namely, p(x) ~ x~ s with 
5 ~ 1.7 (Fig. 10). It is also possible to relate the distribution exponent 5 with dynamic 
exponents (3 and 7 by substituting Eq. (26) into the relation 1 — E(t) = J Xmax dx xn(x, t). 
Using x max (t) ~ the integral is found to behave as t~ 1+f3 ( 2 ~ 5 \ which leads to the 
exponent relation 

(3(2 - 5) = 7. (27) 

Our simulations agree with Eq. (27); e. g., for e = 0, we find (3 = 1, in agreement with the 
obvious intuition that large domains suffer no damage and therefore must grow linearly in 
time. Correspondingly, we find 5 rs 1.7, and 7 0.3. As the fairness parameter e increases 
from to e c , our simulations indicate that the exponents (3 and 7 decrease and both appear 
to go to zero at e = e c . 

The borderline case of e = e c can be expected to lead to marginal behavior which is 
intermediate to the steady state and coarsening regimes. To first establish that e c = 0.5, 
we examined the evolution of a specially prepared system in which a single domain is 500 
times larger than all others. The evolution of this defect domain turns out to be both 
useful and computationally efficient way to ascertain whether the system is in the steady- 
state or coarsening regimes. In the steady state, the defect domain will eventually shrink 
to the average size. Thus e c can be determined as the point where the defect no longer 
shrinks. This approach gives 0.50 < e c < 0.51 with relatively small computational effort, 
suggesting that the value of e c equals 1/2. For e = e c , we do find that the average domain 
age (r(i)) still grows as a power law in time, (r(t)) oc t - 61 , while other basic observables, 
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such as the fraction of empty space E(t), the concentration N(t), and the average size 
(x(t)), exhibit extremely slow variations in time (Fig. 11). On a double logarithmic scale, 
each of these quantities is nearly linear, and a visual fit to the data suggests an exponent 
that is approximately 0.1 or less. A more careful analysis reveals a weak but systematic 
curvature in these data, which suggests that the asymptotic behavior will be slower than a 
power law. However, the time range of the data is insufficient to permit an unambiguous 
fit to a logarithmic or other slowly varying time dependence. 



The existence of power-law domain coarsening for e < e c with a simple relation between 
length and time scales suggests consideration of a deterministic self-similar version of our 
war model to explain the coarsening exponents. We discuss here one such example which 
appears to be particularly suitable for describing the coarsening dynamics for e = 0. In 
this deterministic model, the system starts with a regular array of domains with spacing 
Ax = 1 at t = (Fig. 12). The domain walls move with velocity 1/2 so that the first set 
of warfare events takes place at t = 1. The outcome is defined to be that every second 
domain is annihilated while the remaining domains continue unscathed. This can be viewed 
as arising from an infinitesimal difference in the initial domain sizes. Also at regular time 
intervals At = 1 + , new domains are seeded at the same integer spatial positions of the 
initial domains. Only if the seeding occurs in an empty region does the new domain 
grow. These rules give rise to a pleasing self-similar pattern of domains which resembles a 
Sierpinski gasket, except for the filling of large empty spaces by the continuous input. 

For this system, it is straightforward to compute the properties of domain. These ex- 
hibit strong fluctuations, however, because of the determinism of the model. It is therefore 
convenient to consider quantities which have been averaged over a finite time range, which 
we choose to be between 2 n ~ 1 and 2 n . For example, by low-order enumeration, it is easy 
to verify that between t = and t = 2 n , the total number of domains in existence over a 
length 2 n , starting from the left edge of figure 12, is given by the sequence 2, 7, 22, 67, . . . 
for n = 2, 3, 4, . . .. Solving this sequence, the time integrated density of these domains 
is asymptotically given by |(|) n - If we hypothesize that the density N(t) varies as t~ 7 , 
then the corresponding time integrated density is 



Equating this time integrated density to (2/3) n as determined above, we find the exponent 
7 = 2 — In3/ln2 0.415. In view of the crudeness of this deterministic model, we regard 
this exponent value as being in good agreement with the corresponding numerical result 
of 7 « 0.33. 



(e) Deterministic Self-Similar Model 




(28) 
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IV. SUMMARY AND DISCUSSION 

We have introduced an idealized warfare model in which domains grow at a constant 
rate and where a contact between two domains results in a war, with one or both sides 
suffering casualties. The long time properties of the model are fundamentally governed by 
a fairness parameter e which quantifies the outcome of a war. In a war between domains of 
size % and j > i, the smaller domain is annihilated, while the larger domain emerges with 
a size j — ei. Thus e = 1 corresponds to a fair war where the number of casualties in each 
domain are equal, while e = corresponds to a completely unfair war where the winner 
suffers no casualties. We have examined the long-time kinetics of this model for: (a) the 
heterogeneous process, in an initial distribution of domains is specified, and after which 
evolution by domain growth and intermittent warfare ensues, and (b) the homogeneous 
process, there is a continuous input of infinitesimal domains which then undergo growth 
and warfare. 

For heterogeneous and fair war (e = 1), the system naturally coarsens, with the num- 
ber of domains decreasing as 1/t and their average size growing as t, so that a constant 
asymptotic coverage arises. The domain size distribution obeys scaling in a manner con- 
sistent with these temporal behaviors. While we have not investigated the extension to 
unfair war in detail, the case e = 1 clearly provides a lower bound to the domain size 
distribution for unfair wars. Thus it is evident that the same quantitative linear in time 
coarsening will occur for both fair and unfair wars. 

In the heterogeneous process, there is a wider range of phenomenology which is funda- 
mentally controlled by the fairness parameter. From simulations and a heuristic argument, 
there exists a threshold value e c = 1/2 which separates a steady state regime, for e > e c , 
from a regime of continuous coarsening, for e < e c . In the steady state, the joint age- 
size distribution of domains provides a comprehensive characterization. Interestingly, the 
domain lifetime is relatively large, so that a domain typically survives many wars before 
eventual death. In the coarsening regime, the winner of a war suffers relatively few casu- 
alties which promotes the tendency for the oldest clusters in the system to grow without 
bound. This coarsening evolves according to non-universal e-dependent power laws in time, 
in which the density of clusters and the fraction of empty space decay as t~' y ^\ and the 
average domain size grows as t^ € '. In the completely unfair war case or e = 0, /3(e) = 1, 
as is intuitively clear, but as e — > e c from below /3(e) appears to vanish. At the threshold 
e c , there is a very slow evolution of the system which has yet to be understood. 

The war model also suggests interesting generalizations; a few possibilities and some 
of their attendant consequences are outlined below. 

(i) Size dependent warfare rates: Suppose that the process (z, i+k) — > k occurs at a rate 
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R(i, i + k) which has a power law dependence on the size difference, i. e., R(i, i + k) = k a . 
For fair war, where the losses of each combatant are equal, the rate equations for the 
heterogeneous version of this process are, 

oo oo 

Ck(t) = k^c^+kit) - c k (t) J2\ k ~ ^dit) + Mck-i(t) - c k (t)). (29) 
i=i i=i 

The constant growth suggests the scaling variable x = kj At, but the time-dependent 
prefactor in the domain size distribution may be different than 1/t (compare with Eq. (4)). 
Making the scaling ansatz Ck(t) oc t~ a C(x) and substituting into Eq. (29), self-consistency 
obtains only when a = a + 2. Thus the scaling ansatz is 

c k (t) =t"(«+ 2 ) C(x), with x=^-. (30) 

At 

This scaling form predicts that the size moments behave as power laws in time, M n (t) oc 
t n ~ ol ~ 1 . Hence, the coverage decays with time for positive a, Mi(i) ~ t~ a , indicating that 
aggressiveness which grows with size disparity leads to extinction (perhaps a lesson for real 
civilizations). 

We may further generalize to a size-dependent power-law growth rate A = = . 
In the peaceful limit of growth and no war, the size distribution is peaked around K(t) = 
[(1 + f3)\t] 1 /( 1+/3 \ while warfare produces civilizations in the size range < k < K(t). 
This suggests the scaling ansatz 

CK + /3 + 2 Jc 

c k (t) = t—P^-C(x), with x=— -<1, (31) 

which leads to size moments M n (t) ~ ^( n_a -/ 3_1 )/(/ 3 + 1 ). These results are expected to be 
valid only for (3 > — 1. For (3 = — 1, the typical size grows exponentially, while for (3 < — 1 
the typical size diverges in a finite time, i. e., the most aggressive civilization covers the 
system. 

(ii) Bipolar world: Consider two mutually antagonistic species A and B with aggrega- 
tion occurring when two same-species civilizations (allies) meet, and war occurring when 
dissimilar species meet. If both species grow at the same constant rate, the mean-field 
evolution of A-civilizations is described by the rate equations (in the heterogeneous case) 

oo oo ^ k — 1 

a k(t) = ^6 J (t)a J+fc (t) - a k (t) y^bj(t) + - ^2a i (t)a k - l (t) 

(32) 

oo v ' 

-ak(t) ^2cLj(t) + A(a fc _i(£) - a k (t)), 
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and similarly for S-civilizations. Here a,k(t) and bk(t) are the concentrations of A- and 
S-civilizations of size k at time t. The first two terms on the right-hand side account for 
warfare, the next two terms account for aggregation, and growth is described by the last 
term. 

For this process, it is straightforward to determine that the total number of civiliza- 
tions N(t) ~ 1/t. However, the identification of the appropriate variable in a scaling 
ansatz for the domain size distribution is unclear. In the absence of growth (A = 0), a 
previous study of the resulting aggregation-annihilation process found that the typical size 
scales as k ~ \ft [21]. On the other hand, growth without aggregation and war leads to 
a size distribution which is peaked around k ~ At, while combined constant growth and 
aggregation, but without war [22] , leads to a typical size which grows as t In t . A similar 
ambiguity exists in one dimension since in single-species aggregation and growth, the typ- 
ical size grows as e* [23]. The homogeneous version of the bipolar world model is also of 
interest; a preliminary treatment seems to indicate that a steady state does not arise even 
in the case of fair war. 

(iii) War in two dimensions: It is clearly more realistic to consider our war model in 
two dimensions, where geometrical effects naturally play a more prominent role in defining 
the outcome of a warfare event. If one posits that war is a localized event at the point of 
contact between two domains, then the continued action of war will lead to irregularly- 
shaped domains and possibly to the breakup of countries. These are features which are 
perhaps best investigated by simulations. 

There is a natural simplification which would eliminate the technical difficulties asso- 
ciated with irregularly-shaped domains, however. Namely, start with disk-shaped countries 
and define that after each war the victor retains a disk shape with its center remaining 
fixed. A further simplification is to consider the situation where all domains have the same 
size. In analogy with the corresponding one-dimensional system, analytical progress may 
be possible. For the equal-size case, intuition suggests that the coverage will approach a 
constant value in the long-time limit. Under the further assumption of spatial homo- 
geneity, this would suggest that the number density of civilizations of radius R, N(R), will 
vary as ~ M^/R 2 . 
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APPENDIX A 

We wish to solve the rate equation 



L-^-F(x 7 L) = F(x, L)+x-^-F(x, L) + 6(x-S)G(L) J dyF(y, L)F(x-y-l, L), (Al) 

that describes the dynamics of equal-size growing domains which mutually annihilate when 
they meet. Here F(x,L) = Ln{l,L)/N(L), where n(l,L)dl is the number of neighboring 
civilizations whose centers are separated by a distance which is between [/, l + dl], M{L) = 
f™n(l,L)dl is the total number of surviving civilizations, G(L) = F(x = 1,L), and 
x = l/L. 

To solve Eq. (Al) we apply the Laplace transform, 

/oo 
dxe~ px F(x,L). (A2) 

Note that the relation N(L) = n(/, L) dl can be rewritten as 

dxF(x,L) = l. (A3) 

Combining Eq. (A2) and Eq. (Al) gives 

Eq. (A4) has been solved previously in the scaling limit of L — > oo, where it reduces 
to the ordinary differential equation [14] 

^ = -Cooe- p (l-$ 2 ), (A5) 



whose solution is 



•\>\ P] = tanh | G x / e —dq ) , (AG) 



which contains an as yet undetermined numerical factor Goo. This constant is found from a 
consideration that also establishes the coverage. Civilizations cover the same space, x = 1, 
in units of scaled length, so the coverage is clearly 

Here we use Eq. (A3) and define (x) l by (x) l = dx x F(x, L). 
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In the long-time limit, we use the relation 



oo 



-dq=--\n P - 1 -Y, y -^T, (A8) 



n ■ m 

71=1 



to expand $>(p) in the small p limit as = 1 — 2 exp(2G 0O 7)p 2G °° + — On the 

other hand, from the definition of <E»(p) given in Eq. (A2), we have the expansion, <E>(p) = 
1 —p(x)oo + — Comparing these two forms gives the constant = 1/2 and the (scaled) 
distance between neighboring civilizations, (x)^ = 2exp(7) [14]. This yields the coverage 
in the long-time limit = l/(x) cx> = 0.28073, i. e., Eq. (11). Having established the 
asymptotic coverage, the number density N of these equal-size civilizations asymptotically 
is N(t) ~ Moot -1 in agreement with Eq. (12). 

Now consider the the full time-dependent behavior for which we have to solve the 
nonlinear partial differential equation Eq. (A4). From the form of the asymptotic solution 
of the time independent equation, it is natural to attempt the ansatz 

$(p, L) = tanhftf (p, L)], (A9) 

which allows us to eliminate the nonlinear factor (1 — <E> 2 ). Substituting Eq. (A9) into 
Eq. (A4) gives 



Transforming from the variables (p, L) to u = \pp~L and v = \JpjL simplifies Eq. (AlO) to 

u^-*(u, v) = -G ( -) e~ uv . (All) 
ou \vJ 

The solution to Eq. (All) is now straightforward, 

*(«, v) = e~* v + X (v), (A12) 

up to an arbitrary function x( v )- To determine x( v )i note that the definition of $>(p, L), 
(Eq. (A2)), implies $(p, L) ~ p~ 1 e~ p G(L) in the large-p limit. Since G(L) varies over a 
limited range (as it is clear, e. <?., from relation G(oo) = 1/2), we conclude that $>(p, L) — > 
and hence, *&(p, L) — * as p — > oo. Choose now p ~ L — > oo; in the (it, u) variables, this 
corresponds to u — > oo and i> finite. Thus the integral in Eq. (A12) disappears in this limit 
and we find \l/(oo, v) = x( v ) = implying that x( v ) is trivial. 

Returning now to original variables and replacing £ by rj defined via £ = rju, we rewrite 
Eq. (A12) as 

tf(p,L) = / — G(L v )e- pT1 . (A13) 
A *7 
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We have thus solved Eq. (A2), up to an as yet unknown function G(L). This function can 
be found, in principle, from the initial conditions. Technically, it is convenient to assume 
that there is a finite small-size cutoff L m ; n in the initial distribution which we set to be 
£min = 1 without loss of generality. As an example initial distribution, consider a shifted 
Poisson 

F(x,l) = i e " (x_1) ' x > 1 ^ (A14) 
V ' 1 \0, x < 1. V ; 

In this case 1) = (1 + p)~ 1 e~ p ; therefore, G(L) is determined from the following 
equation: 

drj 



L 



-!- G(rj)e- pT1 = Arctanh[(l +p)- 1 e" p ]. (A15) 
l V 

Although it is impossible to find explicit expression for G(L) in terms of elementary func- 
tions, one can readily compute asymptotics behaviors, e. g., — M(L) oc L -1 . 

Thus the model of equal-size warring civilizations is exactly solvable in one dimen- 
sion. While asymptotic characteristics have been computed by exploiting previously known 
results, the complete solution for arbitrary time is new. However, several related and in- 
teresting properties have not yet been computed. One such quantity is the density of feral 
space, i. e., the fraction of space that has been untouched by any civilization. 

APPENDIX B 

We outline here a solution to Eqs. (20) for the joint age-size domain distribution. 
For this purpose, it proves convenient to rescale length and time by ; c -> x v /2V//j and 
t — > t I V2V fx. In these rescaled units, the previous results for the size distribution become 

n{x) = e —, m Q = N=^. (Bl) 

Then Eq. (20a) simplifies to ^ = -2J(r) + |5(r), whose solution is 

At) = € -^. (B2) 

This shows that in the steady state regime the number density of innocent civilizations, 
I = J °° (irX(r) = |, is one-half of the total number density, / = N/2. Using Eqs. (Bl)- 
(B2), we now reduce the rate equation (20b) to 



[-^ + -^\M{x,t) = 2 dye y M(x + y,r) - M(x,t) 

Introducing g(x,r) = e 3r ~ x M(x, r), reduces Eq. (B3) to 



+ e x ~ 3T . (B3) 



{l + iL) 9 ^ T) = 2 L dy9 ^ T)+i - (B4) 
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Defining now f(x, r) = 2 dy g(y, r) + 1, Eq. (B4) becomes 

One boundary condition, f(x = r, r) = 1, follows directly from the definition of /. To 
obtain a second condition, we compare total number of zero-size civilizations, n(x = 0), 
with the number of innocent zero-size civilizations, X{x = 0) . Both quantities are equal to 
1/2 which means that J °° M(x = 0, t)cIt = 0. Since M(x, r) is nonnegative for all r > 0, 
we conclude that M{x = 0, r) = 0, which leads to the boundary condition f^| x =o = 0. 

Eq. (B5) simplifies further after the change of variables, (x, r) — > (a,/3) = (r, r — x): 

° 2f 2/(a,0). (B6) 



This Klein- Gordon equation is to be solved in the region a > f3 > 0, with the boundary 
conditions 

,1^=1, and m = 0. (B7) 

The symmetry of the governing equation under the exchange of the variables, 
suggests seeking a symmetric solution which depends on a single variable af3. Further 



analysis indicates that the variable z = ^8a/3 is especially convenient. Substituting / = 
f(z) reduces the Klein-Gordon equation to the modified Bessel equation, 

/" + -/'-/ = 0, (B8) 

z 

where the prime denotes differentiation with respect to z. A potential solution is / = 
Io(z). This satisfies the boundary condition Eq. (B6); however, the boundary condition 
of Eq. (B7) is not satisfied. To remove this drawback, we make use of the linearity of the 
governing equation and seek a solution of the form 



f(a,(3)=I Q (z)+(3 n g(z). (B9) 

The first of the boundary conditions in Eq. (B6) is manifestly satisfied (when the index n 
is positive). Substituting Eq. (B9) into Eq. (B6) we get 

g" + 2 ^ ±1 g' -9 = 0. (B10) 

z 

Eq. (B10) is readily solved to find g(z) = Az~ n I n (z), where I n (z) is the modified Bessel 
function of order n. To determine the index n and the amplitude A we substitute Eq. (B9), 
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with g(z) = Az n I n (z), into the second boundary condition in Eq. (B7). This yields the 
relation 

V2h(z) + 2^ J1 AI n _ 1 {z) = 0. (Bll) 
In deriving Eq. (Bll) we used the identities [24] 

2n 

I n -i(z) - I n +i(z) = —I n (z), ^ B12 ^ 

I n _ 1 (z) + I n+1 (z) = 2l' n (z), 

and the equality z = f3\/8 on the diagonal a = (3. Eq. (Bll) shows that n = 2 and 
A = —8. Thus we determine the desired solution to the Klein-Gordon equation with 
mixed boundary conditions 

f{a,P) = I {z)-^-h{z). (B13) 
a 

Returning to the original variables (x, r), and the original joint age-size distribution func- 
tion n(x, r), we get after straightforward computations 

n(x,r) =l(x,r)+M(x,r) = 2xe x ~ 3T ^^ + ^-5(x-t), z = ^8(t-x)t. (B14) 

z 2 

The sum rule Eq. (18a) provides a useful self-consistency check. Substituting Eq. (B14) 
into n(x) = J drn(x,T) gives 

V-V=/ dT W S Ii (Mr - x)t) . (B15) 

1 1 Jx \ 2 (T - X)T V / 



<2{t 

This identity is indeed satisfied [25]. Another quantity which can be calculated exactly 
is average age, (r) = f Q n(r)rdr, where n(r) = JJ" dxn(x,r) (Eq. (18b)). Changing the 
order of integration and using [25], we obtain that (r) = 2. 

Let us finally consider the age distribution of mature civilizations Ai(r). It is given 
by A4(t) = Jq dxA4(x, r). In the large-age limit we use the asymptotic relation [24] 
I\{z) ~ f to estimate the integral. Thus we arrive at 



M(T) - ^ «-<-*>', B = + - 0.977629. (B16) 

Since the age distribution of innocent civilizations decays as e _2r , Eq. (B16) indicates 
that old civilizations are mostly mature. Another interesting computation is the average 
domain size as a function of age, 

, , Sl f n T dxxn(x,r) 

(x{r)) = Jo , T , : (B17) 
J dx n(x, T) 

Using the asymptotic behaviors outlined above, we find 

(x(t)) ~ (2 + 2y/2) (l - . (B18) 
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Table 1 








1 


c 


0.45 


0.46 


0.15 


0.61 


0.4 


0.55 


0.18 


0.59 


0.2 


0.77 


0.25 


0.53 





1 


0.33 


0.35 



Table 1. The estimated values of the exponents that characterize the time dependent properties 
of domains in the coarsening regime, e < 1/2. The estimated error in these numbers 
is 5% or less. 

Table 2 



e 


p(x) 


p(r) 


> 0.5 


exponential 


exponential 


0.5 


1.75 


1.47 


0.45 


1.79 


1.49 


0.4 


1.72 


1.51 


0.2 


1.66 


1.63 





1.67 


1.67 



Table 2. The characteristic exponents of the age and size distributions. The estimated error in 
these numbers is 10% or less. 
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ure 1. Space-time dynamics of vicious civilizations in one dimension. Shown is the hetero- 
geneous version of the war model where the system evolves from a fixed initial state. 
When two civilizations of size i and j > i meet, the result is a diminished civilization 
of size j — ei. By definition, casualties occur at the point of contact, so that the side 
of the large civilization on the battle front retreats by e i while the other side is not 
affected. 
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Figure 2. Space-time evolution of two domains for the special case of e = 1/2. The position 
of the frontier of the winner immediately after the war and at the birth-time of the 
smaller combatant is the same. 
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Figure 3. Scaled civilization size distribution, n(x,t)/N(t), from numerical simulations of 500 
configurations of heterogeneous war with 10 5 initial domains in one dimension. The 
data are for t « 1.5 10 « 57.7 (o), t « 1.5 12 « 129.7 (+), and t « 1.5 14 « 291.9 (*). 
Also shown is the corresponding mean-field result ^ sin ^| (dashed). 
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Figure 4. Schematic illustration of the equivalence between the space-time evolution in the one- 
dimensional war model with equal size domains, and the "dual" problem of coarsening 
of domains through the successive elimination of the smallest domain. In this equiv- 
alence, the centers of the domains in the war model are identified with the positions 
of domain walls in the coarsening process (dashed lines). The two domains which are 
eliminated in the dual problem are indicated by the dark shading. 
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Figure 5. Correspondence between various interaction events and the terms in the rate equation, 
Eqs. (14). (a) Illustration of the first term on the right-hand side of Eq. (14a), (b) the 
third and (c) the fourth terms on the right-hand side of Eq. (14b). 
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Figure 6. Steady state domain size (o) and age (□) distributions for the case e = 1. The corre- 
sponding predictions from the solution to the rate equations are also shown (dashed 
lines) . 
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Figure 7. The r-dependence of the steady-state size (*) and age (o) correlation functions, C s (r) = 
(siSi +r ) / (si) 2 — 1 and C T (r) = (TiTi +r ) / (n) 2 — 1, respectively, for the case of fair war 
e = 1. vfill 
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Figure 8. Average domain size as a function of domain age in the steady state for the fair war 
case of e = 1. Shown is are the simulation results (□) and the predictions based on the 
solution to the rate equations (solid line). 
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Figure 9. Representative simulation results for time-dependent quantities for a system that has 
reached a steady state, e = 2/3 (o), and for a system which perpetually coarsens, 
e = 2/5 (+). Shown are: (a) the total number of domains, (b) the average domain 
size, (c) the fraction of empty space, and (d) the average domain age. This data is 
based on 20 configurations of system with 10 3 domains initially 
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Figure 10. Domain size distribution in the coarsening regime for e = (o), e = 0.2 (+), e = 0.4 
(•), and e = 1/2 A). To distinguish the different data sets, the points for each value 
of e have been shifted vertically by a small fixed amount. 
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Figure 11. Simulation results for time-dependent quantities for the marginal case of e = 1/2 
Shown are: (a) the total number of domains, (b) the average domain size, (c) the 
fraction of empty space, and (d) the average domain age. This data is based on 100 
configurations of system with 10 3 domains initially. 



41 




Figure 11 (d) 
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Figure 12. Illustration of the deterministic version of war dynamics for the case e = 0. Domains 
are created at regular lattice points and at regular unit time intervals. When equal 
size domains meet, a deterministic rule is implement to decide the survivor. 
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